MIT-CTP-2501 
hep-ph/951243S 



Supernatural Inflation: Inflation from 
Super symmetry with No (Very) Small 

Parameters 

Lisa Randall Q 
Marin Soljacic [] 
Alan H. Guth fl 

Laboratory for Nuclear Science and Department of Physics 
Massachusetts Institute of Technology 
Cambridge, MA 02139 

Abstract 

Most models of inflation have small parameters, either to guarantee sufficient inflation 
or the correct magnitude of the density perturbations. In this paper we show that, in 
super symmetric theories with weak scale supersymmetry breaking, one can construct viable 
inflationary models in which the requisite parameters appear naturally in the form of the 
ratio of mass scales that are already present in the theory. Successful inflationary models can 
be constructed from the flat-direction fields of a renormalizable supersymmetric potential, 
and such models can be realized even in the context of a simple GUT extension of the MSSM. 
We evade naive "naturalness" arguments by allowing for more than one field to be relevant 
to inflation, as in "hybrid inflation" models, and we argue that this is the most natural 
possibility if inflaton fields are to be associated with flat direction fields of a supersymmetric 
theory. Such models predict a very low Hubble constant during inflation, of order 10 3 -10 4 
GeV, a scalar density perturbation index n which is very close to or greater than unity, and 
negligible tensor perturbations. In addition, these models lead to a large spike in the density 
perturbation spectrum at short wavelengths. 
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I. INTRODUCTION 



Inflationary models M in general require small parameters in the particle theory La- 
grangian, to provide the flat potential needed for sufficient inflation and for the correct 
magnitude of density fluctuations. The need for unmotivated small parameters tends to 
weaken the credibility of a theory, so one hopes that the origin of these parameters can be 
understood. It is conceivable, of course, that the explanation lies beyond our present under- 
standing, just as we presently have no accepted explanation of why the Yukawa coupling of 
the electron is 2 x 1CU 6 , or why the weak scale lies 17 orders of magnitude below the Planck 
scale. Nonetheless, it would be encouraging to find that the small parameters required by 
inflation could be obtained from small parameters that are already essential to the particle 
theory, so that no additional small parameters are introduced. Models in which the small 
parameters arise as ratios of known particle physics mass scales are particularly attractive 

From a field theoretical perspective, it is difficult to see how flat direction fields can be 
present in a nonsupersymmetric theory (with the exception of Goldstone bosons, considered 
in Ref. ||) given the effect of radiative corrections . We will assume therefore that the 
world is super symmetric, and ask whether an inflationary potential can arise naturally in 
the context of the mass scales which we expect might be present. Some natural candidates 
for these scales could be the Planck scale M p « 10 19 GeV, the GUT scale Mq « 10 16 GeV, 
the intermediate scale Mj 10 11 GeV, and the supersymmetry breaking scale 7713/2 ~ 1 
TeV. 

However, straightforward considerations show that it is difficult to implement this strat- 
egy. The magnitude of density perturbations points to the GUT scale as setting the energy 
density during inflation, since (M G u T /M p ) 2 ~ 10~ 6 . Although suggestive, it is difficult to 
exploit this high scale in an inflationary model. In Ref. 0] it was argued that the necessity 
for cancelling the cosmological constant after inflation provides significant restrictions on a 
model in which the inflation scale is greater than the supersymmetry breaking scale. Models 
have been suggested in which the energy density in the early universe is well above the low 
energy supersymmetry breaking scale, and might be governed for example by the value of 
a moduli field 0H. For these models, it must be checked that the constraint of Ref. || 
is satisfied. Beyond this, however, it is difficult to study such models in detail without a 
concrete realization. 

It might also be that the inflation scale is generated by GUT physics. An interesting 
example of this type of model is Ref. |7|] . It is however questionable whether the GUT scale 
exists as a fundamental scale of particle physics at all in light of the doublet-triplet splitting 
problem J8|. Furthermore, it is likely that the inflaton field would be charged under the GUT 
group and the reheat temperature would be too high (in excess of the gravitino bound). 

In the context of supersymmetric models, an attractive scale for the vacuum energy 
density during inflation would be set by the intermediate scale, Mj. This scale is very 
likely to be present in a hidden sector model of supersymmetry breaking. It is also the 
right energy scale for the potential associated with moduli fields, which might be natural 
candidates for flat directions. The problem here, however, is that simple dimensional analysis 
arguments (to be reviewed in Sec. 2) show that density perturbations would generically be 
either far too small or far too large, depending on assumptions. For the case of inflation 
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driven by a single moduli field, the dimensional arguments show that the requirements 
of sufficient inflation and correct density perturbations imply that 1) the variation of the 
inflaton (moduli) field during inflation is of order M p , and 2) the energy density during 
inflation is of order M 4 where M ~ 10 16 GeV. An energy density of order Mf would 
produce density perturbations too small by about ten orders of magnitude. For the case of 
a chaotic inflationary scenario, the variation of the inflaton is again of order M p , but in this 
case the density perturbations are much too large unless the quartic coupling A is about 

io- 12 . 

In this paper, we show that the argument that inflation at an intermediate scale is 
untenable lacks sufficient generality, and can evaporate if one drops the assumption that 
inflation is driven by a single scalar field. We describe a class of two-field models, for which 
dimensional analysis estimates show that 1) the variation of the inflaton is of order Mi or less, 
and 2) the energy density during inflation is of order Mf. We then go on to illustrate these 
ideas with models motivated by supersymmetry with soft supersymmetry breaking. We will 
find that these models not only solve the naturalness problem of obtaining sufficiently many 
e-foldings of inflation, but also generate very nearly the correct size of density perturbations 
based on the parameters of supersymmetry breaking. We therefore refer to our models as 
"supernatural" inflation. 

This model contains a similar structure to the "hybrid" inflation models, proposed by 
Linde and studied by Copeland, Liddle, Lyth, Stewart, and Wands [[UJ. The fact that the 
standard dimensional naturalness arguments for the number of e-foldings and for Sp/p do 
not apply, and that the Hubble scale during inflation will be low was also clearly recognized 
by these authors. Our point here is to emphasize that the most natural scales for successful 
implementation of two field inflation of the "waterfall" type are the scales associated with 
supersymmetry breaking and the Planck scale. Furthermore, our models more accurately 
reflect masses and couplings associated with flat direction fields, and we will motivate the 
parameters and potential we use by consideration of flat directions in the MSSM. Hybrid 
inflation in the context of SUSY leads one to the interesting conclusion that the Hubble 
scale during the inflation which established the density perturbations might have been of 
order 10 3 -10 4 GeV, rather than 10 13 GeV. 

In the following section, we present the general arguments for why supersymmetry scales 
do not work in single field inflation models. We then review the general idea of "hybrid" 
or "waterfall" [|l(J models, and show why the single-field arguments do not apply to the 



two-field case. In Section 3, we present supernatural inflation models, in which we assume 
the inflation sector consists of flat direction fields whose potential is generated through 
supersymmetry breaking and nonrenormalizable operators. We derive the constraints on 
parameters consistent with the requisite number of e-foldings and density perturbations. In 
Section 4 we explore the possibility of a renormalizable coupling between the flat direction 
fields. In the following section we motivate the models of Sections 3 and 4 by briefly con- 
sidering flat direction fields in the standard model, and present an illustration of the model 
of Section 4 in the context of a GUT extension of the MSSM. In Section 6, we give details 
of the evolution of the two fields in our models. We analyze the density perturbations that 
result from this evolution, and discover a novel spike that is predicted to appear at short 
wavelengths. Such a spike could lead to overproduction of black holes, but we show in Sec- 
tion 7 that existing constraints on black holes are satisfied for the parameters of interest. In 
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Section 8, we show that the thermal production of gravitinos is also not a problem in our 
models. In the following section we show that successful baryogenesis can be accomplished 
in the context of late inflation. In Section 9, we discuss miscellaneous aspects of our models, 
and in Section 10, we conclude. 



II. ONE VS. TWO FIELD INFLATION 

We begin this section by reviewing the "standard" arguments for why the inflaton in 
"natural" inflationary models varies on the scale M p and why the scale for the energy density 
should be larger than the intermediate scale in inflationary models with a single field. 

For the purposes of these dimensional arguments, we first assume the potential takes the 
form 

v = M'g^/f) (i) 

where Q is a bounded function of order unity. Here we have in mind for example a moduli 
field, with M ~ Mj. If we assume the slow roll equation of motion 3H<f) = —V, where H is 
the Hubble constant during inflation, the number of e-foldings is 

N= [Hdt= = (2) 

J J (f> J MjV M p M p Q> v ; 

There are essentially two possibilities. If Q is a bounded function, and Q' is not very tuned 
to have very flat sections, one is in the regime of what might be expected for a moduli type 
field. In this case, the requirement of about 60 e-foldings of inflation favors / of order M p 
and a change in during inflation at least of order M p . Even when this is satisfied, some 
tuning of the potential is required. 

The alternative possibility is that one is in a chaotic JTT[ inflationary scenario, in which 



case Q will be dominated by monomial behavior for sufficiently large field, and V ~ V/0. 
In this case, / is not defined, but one would still conclude A0 rs M p . 

Density fluctuations are also readily estimated under the assumed form of the potential. 
They are given by 

p ~ ~ V ~ \M P ) M p Q> y } 

Assuming a potential of the moduli type, with Q and Q' of order unity and / of order M p , 
we find that ^ ~ (rf^J favoring M ~ 10 _3 M p . Detailed calculations might change M by 
an order of magnitude or so, but it is clear that M ~ 10 11 GeV ~ Mj is strongly disfavored. 

In a chaotic scenario on the other hand, one would conclude that the density fluctuations 
are too large unless there is a small parameter. For example, a simple dimensional argument 
would lead to the conclusion that for V = \<p A , A ~ 10~ 12 . Without further motivation for 
these small numbers, such a potential seems unlikely. 

So one is led to the conclusion that it is difficult to naturally obtain sufficiently many 
e-foldings and the correct magnitude of density perturbations, without invoking either small 
numbers or a new mass scale. 
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It is apparent, however, that there is a loophole in the above argument. From Eq. (0) 
it is clear that the constraints on / and A0 during inflation arise because it is the same 
potential V(0) that controls the inflation rate H and the speed of the inflaton field 0. These 
constraints can be avoided, therefore, if the energy density during inflation is provided from 
some source other than the scalar field which rolls and controls the ending of inflation. 

The simplest way to implement this idea would be with two fields. This idea is essentially 
that first proposed by Linde |IIJ as "hybrid" inflation or "waterfall" models. There are two 
fields if) and 0. The first field, which we call the inflaton, has a very flat potential. It starts 
at a large field value, and slowly rolls (via its classical field equations) to the origin. 

The second field, 0, has a potential whose minimum is far from the origin. In most 
previous incarnations of hybrid inflation, the scale of variation of this field is Mi, though 
in our models the scale will be M p . When if) has large field value, it gives a positive mass 
squared term in the potential at the origin, so the classical field equations push to the 
origin. When if) gets sufficiently small (of order Mj or less in our models), the mass squared 
of goes negative, and makes the transition from the origin to M p . 

The key feature of this model is that the energy density during inflation is dominated by 
the potential energy of the field at = 0. There are no tunings in the if) potential to get 
a small mass during inflation and a large mass afterwards, since its mass is always small, as 
is its potential energy. Because H depends on the value of and is not determined by the 
field if), which acts as a switch to end inflation, the naive estimates do not apply. 

The second key feature of this model is that the ending of inflation is controlled by when 
the mass squared at the origin changes sign. One can obtain a large number of e-foldings 
with the variation of the inflaton field ip much less than M p . Let us see this explicitly. 

We assume a potential which takes the form 

V = M 4 £(|0|//)+s(|0|,H) + m 2 H 2 (4) 

where M and / are to be determined, the function g is the term responsible for the if) 
dependence of the mass, and m is of order 777.3/2- 
We now have 

J if) J m 2 lf) m 2 M 2 \lf) fi na i) ^ ' 

Notice that the scale M in the numerator is independent of the mass and coupling of the if) 
field (in the limit that the if) contribution to the energy density is small) so that the previous 
arguments for one-field inflation no longer apply. Clearly for inflation to give several e- 
foldings requires only that ip changes by an order of magnitude, and that M 4 > m 2 M 2 . No 
inflaton variation of order M p is required, and so far, it seems M ~ Mj could be a good 
choice. 

Let us now consider density fluctuations under the same assumed form for the potential. 
We find 

5p H 2 H 3 M 6 

JL ~ ~ ~ (6) 

p ip m 2 if) M 3 m 2 ip 

The point is that the numerator H 3 has its scale set by the potential energy while the 
denominator is determined by the if) field. We construct a model so that if) at the end of 
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inflation is of the order Mi or smaller. If we also take M ~ Mi, we find Sp/p ~ Mi/M p 
or bigger (rather than (Mi/M p ) 2 as was the case in single field models). Although the 
coupling between the fields can have a coefficient which varies by many orders of magnitude, 
as does if) in eqn. (|6|), the strong M dependence of Eqn. (§) allows for agreement with 
the COBE constraint with only a relatively small M variation. This is very promising from 
the perspective of relating inflation models to real scales of particle physics. To answer the 
questions of how well these ideas really work, and how constrained the parameters of the 
models really are, requires a detailed investigation of particular examples of these ideas. 



III. SUPERNATURAL INFLATION 

We define Flat Direction Hybrid Inflation (FDHI) models as those motivated by the 
properties of moduli fields or flat directions of the standard model. For moduli fields with 
no gauge charge or superpotential, the whole potential arises from the Kahler potential once 
supersymmetry is broken. This potential for and if> wm take the form 

V(<P^)=Mff{<P/M p ^/M p ) (7) 

where the dimensionless coefficients in / should be of order unity. 

However, it is clear that a model of this sort will not give rise to inflation with sufficiently 
large density fluctuations, since during the relevant period if> will typically be of order M p , 
and the resulting Sp/p will be of order (M//M p ) 2 . We conclude that it is essential to have 
an additional interaction between if) and 0. In this model, we assume the existence of a 
superpotential which couples if) and but which is suppressed by a large mass scale M' . For 
standard model flat directions, such higher dimension operators are to be expected, with 
M' equal to M p , Mq, or some dynamical scale. In the case of moduli fields, it might be that 
this scale is of dynamical origin; one can readily determine how the answer changes with 
the form of the superpotential and the size of the mass scale. 

We therefore assume the presence of a superpotential. The example we take is 

6 2 ib 2 

W = z-f- (8) 
2M' K ' 

We now need to specify the form of the supersymmetry breaking potential. We assume 
both ip and have mass of order the soft SUSY breaking scale of order 1 TeV (where we 
will need to test the consistency of this assumption). We assume that the potential for 
the if) field gives a positive mass squared at the origin, while the field has negative mass 
squared at the origin. Furthermore, we assume that the cosmological constant is zero at the 
minimum of both if; and 0. The specific form of the potential we choose is 



, / l<&l\ n n l$| 4 |\l/| 2 + l\P| 4 |$l 

V = MW ( J|ij + ml |vlf + m 1 <0> 



where again we assume (and verify for consistent inflation) M ~ Mj. When the parameters 
are motivated by supersymmetry breaking, we refer to our models by the name supernatural 
inflation. We will see that one very naturally obtains the correct magnitude of density 
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perturbations, and sufficiently many e-foldings of inflation, using parameters and a potential 
which are well motivated in supersymmetric models. 

For the purpose of the inflationary model, the scalar field can be assumed to be real. 
Defining \I/ = (if) + iipi) / V% and a similar equation for $, the potential for the real fields 
becomes 

^W(^2/) + ^ + ^W (10) 

The if) mass is and the magnitude of the (imaginary) mass term (at the origin) is 
m$ = M 2 /f. During inflation, (f> is confined near the origin. The field if) slowly rolls towards 



the origin and inflation ends about when if) = ip c = ^j2M'm$. It will turn out that either 
fn^j/m^ or M'/M p is small, so that during inflation the term m^if) 2 is small relative to M 4 . 

The Hubble parameter during inflation is therefore approximately H = y87r/3M 2 /M p . 

We expect / is of order M P , or equivalently, m| is of order m^ 2 . Although it looks 
like we took a very special form for the (^-potential in Eq. (Q), the use of the cosine is not 
essential. As can be seen from a Taylor expansion, only at the very late stages of inflation 
are terms other than the constant and mass term relevant. We could equally well have 
specified a potential which is truncated at fifth order in the fields, or which has different 
higher order terms. Although both if) and <f> might be moduli or standard model flat direction 
fields, we assume their potentials are of very different form; the particular case we assume 
is illustrative of how a model could work. 



The constraint from density perturbations in the slow-roll regime is [|12] , p~3 

y3/2 



M*(dV/diP) 

where M p = M p f \/%tx. This gives the constraint 



6-10" 4 (11) 



where 



_4^-W^e riV = 6.7-10- 6 (12) 
mlMl V M' y 1 



' —l + Vf + A^lf (13) 

where the approximation in Eq. (|13|) is required if the slow roll conditions are satisfied. 
Here we have defined ^ = m^/H and have measured time in e-foldings away from the 
time N = when if) = if> c (where inflation ends at positive N). It is clear that a lower 
M' makes the value of if) at the end of inflation lower, which in turn increases the density 
perturbations. The exponential in Eq. ( |1"2"D determines the scale dependence of the density 
perturbations, characterized by the scalar index. 

The scalar index a s is readily determined from the scale dependence of the density per- 
turbations to be — /i^/3. This can be seen directly from the formula for density perturbations 
above. Alternatively, it is extracted from the general formula [13] 
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n = 1 — 2a 



s 




(14) 



where dimensionful factors should be compensated by M p . Notice that the second term is 
negligible for all models for which the inflaton field value is much less than M p . This is 
readily seen from the fourth expression in Eqn. which implies M p V'/V ~ A<f> / (N M p ) . 
The third term is positive in our model, because the inflaton field rolls toward, rather than 
away, from the origin during the end of inflation. 

We see for this model that n is always greater than 1, and is very close to 1 for small 
fi^p, which is the case for large M'. This differs from the usual prediction for new inflation 
or chaotic inflation models. The current upper bound on n is uncertain as is summarized in 
Ref. jn|. These bounds, along with the validity of slow roll, prevent too large values of 
From Figures 1-3, we see large fi^ is permitted only for the smallest value of M', where the 
bound on n will provide an additional constraint. 

The fact that n is greater than or close to 1 is a characteristic feature of our models, 
which should help them to be distinguishable in the future, when a good measurement of 
the CMBR is obtained. 

Another distinctive feature of these models is that the ratio of the tensor to scalar 
contribution to the quadrupole 



Again this follows from the small value of the inflaton field ip near the end of inflation. 

As we have argued in the first section, models of inflation which have only a single field 
should have the inflaton field taking a value of order M p near the end of inflation if 50 
e-foldings are to be obtained without fine tuning. The combination of negligible 1Z and n 
never below 1 are distinctive features of these models which should help distinguish them 
from other possible inflationary models in the future. 

In Figs. 1-3, we show values of the parameters when M' ~ M p , M' ~ Mqut and 
M' m Mi respectively. The values shown were found by imposing the correct magnitude of 
density fluctuations and choosing the minimum = m^/H consistent with a sufficiently 
rapid end of inflation (see Sec. 6). We chose the range of M to optimize parameters. Smaller 
M would increase the values of l/fty and fi^. Large M would improve (that is, decrease) 
these ratios but would make the masses uncomfortably large relative to the TeV scale. We 
find that smaller M' gives more natural ratios for the mass to the Hubble scale, though in 
all ratio of less than 100 can be obtained. 

These constraints assumed that the contribution of <fi to density perturbations was small. 
In order to check the consistency of this assumption, we need to consider the evolution of 
if) and in the late stages of inflation. It will turn out that inflation must end reasonably 
quickly after if> reaches ip c so that perturbations exit the horizon while the field is still 
confined to the origin. This gives a lower bound on /i^,. In Sec. 6, we will investigate the 
mass constraint in detail. 




(15) 
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FIG. 1. Parameter choices for supernatural inflation with M' at the Planck scale. Here fi^p is 
chosen to give the correct magnitude of density fluctuations for the minimum fi^ consistent with 
a sufficiently rapid end to inflaton. 
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FIG. 2. Parameters for supernatural inflation with M' at the GUT scale. The values were 
chosen by the same criteria used in Fig. I. 
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FIG. 3. Parameters for supernatural inflation with M' at the intermediate scale. 



IV. ANOTHER MODEL 



In Section 3, we investigated the possibility that there is a nonrenormalizable superpo- 
tential. However, it is frequently the case that flat directions lift each other; that is, the 
renormalizable potential does not permit certain field directions to be simultaneously flat. 
In this section, we present an alternative model with a renormalizable potential. It will 
turn out that this model requires a small coupling in the potential. We will motivate this 
assumption in Section 5, where we consider particular choices of ip and <fi chosen from the 
supersymmetric standard model where we will show that the small coupling can actually be 
related to one of the known small Yukawa couplings! Once we assume this small parameter 
(again an unexplained but perhaps necessary parameter of the MSSM) we will find that /i^ 
and can both be close to unity. 

So we take the potential to contain the soft supersymmetry breaking terms as before but 
to contain a renormalizable coupling between ip and <fi. Specifically 



This model has the essential features of the FDHI model of the previous section. The 
difference is the value of ip c which in this model is 



The density fluctuations give the constraint 




(16) 
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FIG. 4. Parameter choices for supernatural inflation with renormalizable couplings, for 
A = 10" 4 . Like the previous figures, this graph shows the parameters associated with the minimal 
allowed value of 



where 




9,2 f4 



(19) 



If A is of order unity, to satisfy Eq. (|T8j) requires /i^ > 10 . However, if A ~ 10 — 
10 -5 , the model works perfectly with /i^, and fi^ both of order unity. In Fig. 4 we plot 
the parameters of the model for A = 10~ 4 , where again we have chosen the minimum fi^ 
consistent with a sufficiently quick end to inflation. We see there is virtually no fine tuning, 
so long as a small A exists(!). In the following section, we explain why such a small value 
for A is not necessarily unexpected. 



V. EXAMPLES OF "FLAT DIRECTION" INTERACTION POTENTIALS 

Up to this point we have considered abstractly how to make a successful inflation model 
premised on properties motivated by flat direction fields. In this section, we motivate the 
sort of models we have considered by demonstrating examples of flat direction fields in the 
MSSM whose couplings are those required for a successful supernatural inflation model. 
However, because the field minimum is not at zero, which in the context of the standard 
model would imply large gauge symmetry breaking, the minimal standard model is not 
appropriate. Nevertheless, in simple extensions of the standard model, there are neutral 
fields which can play the role of the <ft field. We will present an example of a GUT extension 
of the MSSM which could contain appropriate "i/j" and "0" fields. 

We first present a few interesting cases of the form of the superpotential which arises in 
the context of the MSSM. Other flat direction possibilities could motivate further general- 
izations of the models we have studied. 
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Flat directions of the MSSM were considered in Ref. [ff6] and a complete analysis was 



given in Ref. ||17|| . They can be parameterized by gauge invariant combinations of fields. 
For example, the flat direction H\ = ifj = will be denoted as <fi = H\H\ where we have 
given explicit gauge indices. In each case, it is necessary to explicitly check for F-flatness 
by examining the form of the superpotential. 

Example: ip = ui l d\d\, <p = Qld^L 2 W = QuQd/M' . The superscripts here refer to 
the color indices and the subscripts to flavor indices. We haven't specified the indices in 
the superpotential where in principle all allowed contractions can appear in different terms. 
This model would agree with the physics of our first model of inflation, except that the 
superpotential actually takes the form tp 2 <fiX / M' + a<j) 2 ipX/M' , where X is a field which is 
not flat and which stays at zero during the relevant stages of inflation. This can be seen 
by explicitly substituting u\ = d\ = d\ = ip and Q\ = d\ = L 2 = 4> into all possible 
terms of the above form in the superpotential (that is with arbitrary flavor and allowed 
gauge contractions). The potential will be as we have studied, except that the coefficients 
of the ip 4 (f) 2 and 4 ?/> 2 terms need not be identical, since they arise from distinct terms in the 
superpotential, and there is a cross term which changes the exact evolution of <fi but in no 
significant way, which arises because the separate superpotential terms can depend on the 
same X. 

Example: ip = LH U , <ft = H u Hd 

In this example, SU(2) D terms will lift the flat direction. This can be seen by solving 
for the fields in terms of the flat directions and substituting into the D terms of SU(2). It 
can be seen that the D term does not vanish, but can involve a ip 2 <j) 2 cross term, suppressed 
only by g 2 , where g% is the SU(2) gauge coupling. As discussed in Section 4, a model of this 
sort with a large gauge coupling can work, but requires tuning 

Example: ip — udd, (ft = LH U , W = X u Q u H u u 

This example realizes perfectly our scenario with a renormalizable potential generated 
by a small Yukawa coupling, but not a large gauge coupling. In fact, if it is indeed the up 
quark in the ip field, the potential works as well as could be hoped, since it depends on 
X 2 ^ 2 ^ 2 which as we have shown gives fi^ and ^ of order unity. In this model, the small 
size of density fluctuations arises as a natural consequence of the small up quark Yukawa 
coupling! 

There are in fact other interactions in this model, with the field Q u . However Q u is 
not a flat direction and is assumed to be zero (or small) throughout inflation, so that it is 
irrelevant to the analysis. 

In fact there are many examples of the above type. Even with somewhat bigger Yukawa 
coupling, the correct magnitude of density fluctuations can be obtained at the expense of 
a larger ratio of fi^. This is probably the nicest possibility for realizing the inflationary 
scenario we have outlined, because the only small numbers are those already present in 
the form of Yukawa couplings. There are no unlikely assumptions required for the correct 
magnitude of density perturbations and a sufficiently rapid exit to inflation. 

The problem with the MSSM as the source of inflaton candidates, as we have already 
stated, is that <fi has a nonzero expectation value at the end of inflation. Because 4> in 
general carries standard model gauge charge, this is not permitted. However, in GUT or 
other generalizations of the MSSM (or in models with completely independent fields which 
do not carry standard model gauge charges) one can readily realize the scenario we outlined. 
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Model: Consider a generalization of the standard SU(5) GUT theory to SU(6), where 
now the Higgs fields are in the 6, 6, and 35 (H, H, S) representations of SU(6). The matter 
consists of three generations of 15 + 6 + 6'. This model has been well studied in the context 
of solving the doublet triplet splitting problem |18|]. Examples of specific models with the 
requisite accidental symmetry were presented in Ref. [[H] . 

Here we assume there is no superpotential for the H and H fields, but that the potential 
created by the soft supersymmetry breaking terms is minimized at (H) ~ (H) ~ M p . 
Notice that this will break SU(6) to SU(5) which can survive to the GUT unification scale, 
and is therefore phenomenologically consistent. The £ field acquires a vacuum expectation 
value of order Mq. We assume that the £ field acquires this expectation value through 
renormalizable interactions, and is therefore not flat, and furthermore has reached its true 
minimum at the time of inflation. 

Now consider = H°H° and ip — 152 2 6f 6'g where we have labeled the matter according 
to its generation number (subscript) and according to the SU(6) index (superscript) (where 
indicates the SU(5) neutral direction). 

The superpotential which is required to complete the model can readily be chosen in 
accordance with the requirement of a small Yukawa coupling. To explicitly write the term 
is subtle however for two reasons. First, the leading (renormalizable terms) which give mass 
are the terms which make the 6' fields heavy and the term which gives the top quark a mass. 
However to give a renormalizable top quark coupling requires that the top be in a 20 of 
SU(6). All other masses arise from nonrenormalizable operators, and therefore appear more 
complicated. However, because of the large expectation values of the S, H, and H fields, 
these terms reduce to ordinary Yukawa couplings. 

A toy model which would give the necessary Yukawa coupling would be W = A A1 15 2 -ff6 2 , 
where A M ~ 1CT 3 . This example resembles the up quark example. However, this is only a 
toy model because such a term actually gives the 6' a mass, and in fact defines the 6' fields. 
If this Yukawa coupling happens to be small, the density fluctuations would be small. Since 
we know little about the extra quark Yukawa coupling, we give a model involving the known 
quark mass parameters. 

The higher dimension term W = V& c (20£)(-^)152 can generate the mixing angle between 
the second and third generation. The effective Yukawa coupling between the flat direction 
fields from this term is Vb c E/M p which is about 10~ 4 (since the £ VEV breaks SU(5) to 
SU(3) x SU(2) x U(l) at the GUT scale). The density fluctuations in this model are then 
naturally of order 10~ 4 . 

The last model works very well, as is illustrated in Figure 4. It is extremely interesting 
that the inflation scenario we have devised can be explicitly realized in the context of a 
known model of particle physics. 



VI. AND EVOLUTION 

In this section we will discuss the details of the evolution of ip and 0. We can then derive 
the constraint on the mass. 

As we have argued, inflation ends at about the time the squared mass changes sign at 
the origin so that will roll towards its true minimum. However, because the mass is not 
large compared to H as in previous implementations of the hybrid inflation scenario ||10|| , a 
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careful study is required to ensure that inflation ends sufficiently rapidly that our formula 
for density perturbations applies. We will see that if inflation ends too slowly, perturbations 
will leave the Hubble radius when the fluctuations are large. In this scenario, if these 
fluctuations were on measurable scales (say greater than 1 Mpc and smaller than 10 4 Mpc) 
either the size of density fluctuations or the deviation from a scale invariant spectrum would 
exceed the experimental bound. By a detailed study of the combined evolution of the -0 
and fields, we determine the necessary constraint on the mass for consistency of our 
model. However, throughout this section, it should be remembered that this constraint is 
only very important when /i^ is small, because it slows the transition which causes the end 
of inflation. When /i^ is close to unity, values of /i^ near unity are also adequate for a rapid 
end to inflation. For this reason, we will focus in the discussion here on the case of small /i^,. 
We present in detail the analysis for the model of Section 3, where the constraint is more 
severe. A similar analysis was done for the model of Section 4 in order to obtain Figure 4. 
We first consider the time evolution of -0. The equation of motion for -0 is 

^ + 3^ + m 2^ + J^ = (20) 

where the dot denotes derivative with respect to time t. 

There are three relevant stages of evolution of the -0 field. In the early stage of its 
evolution when the field is small (and so is 4 /M' 2 ), the -0 field obeys the slow roll 
equation of motion as it evolves towards the origin. 

i/>(t) = ij c e~^ N/3 (21) 

where ip c = ^jlM'm^ is the value of -0 when m^(N = 0) = (at the origin) where we 
measure time in units of if -1 . Eventually, the field will grow to a sufficiently large value 
0c where the -0 mass becomes large, and the -0 field acts as a coherent state of oscillating 
particles with mass m$(t) = 2 /2M'. (Here we use the argument t to distinguish the time- 
dependent physical mass of -0 particles from the time-independent mass parameter in the 
potential.) We define C by <p 2 c /2M' = H. 

In our numerical simulations, we replace the time evolution of -0 by the time evolution 
of its envelope at a time sufficiently late that we can neglect the term -0 4 2 /8M' 2 in the 
potential. The envelope ip e obeys the approximate equation of motion 

4 = - + \m) ^ (22) 

where T(t) is the -0 decay rate, where the time dependence arises from the time dependent 
mass. Without further knowledge of the identity of the -0 field the decay rate is an unknown 
parameter of the theory. We constrain the model under two reasonable scenarios for the -0 
decay. If -0 has renormalizable couplings to other fields, its decay rate can be as large as 
Tb ~ m^(t). Of course there are unknown coefficients to this estimate but this probably 
represents the maximum possible rate. The true rate should lie between T b and T[, where 
Ti = m^(t)/(M p /v / 87r) 2 . This latter decay rate assumes no renormalizable couplings, but 
Planck suppressed interactions which allow the -0 field to decay. We evaluate the final stages 
of evolution allowing for these two possibilities for the decay rate. After the time at which 
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~ H, the amplitude of the if) envelope is quickly reduced to zero, and it is only the (f> 
field which remains. Notice that m(t) could exceed M p before the end of inflation. However, 
the if) field decays well before this inconsistency in the expansion is reached. 

Because of the strong dependence on the if) field of m^it) 2 = m|(l — ip 4 /ip^), the if) 
evolution is critical to determining the <f> evolution. In the early stage of the <f> evolution, 
it can be described by a Fokker Planck probability distribution P(<f),t). This distribution 
will be centered at the origin, but the spread will determine the effective amplitude of 
the field 6 = 



Eventually the effective amplitude will be sufficiently large that the 
classical equations of motion take over. By considering the exact solution to the Fokker 
Planck equation, we determine the correct initial condition for the subsequent evolution of 
the classical field equations. 

The Fokker Planck equation in a de Sitter background with time independent Hubble 
constant H is [EDI 



dP((f),t) H 3 



dt 



d 2 



m 2 (t)(f) 
3H 



(23) 



where we have assumed the slow roll equation of motion to be valid. With the evolution of 
if> described by Eq. (E_Tj), the time-dependent mass of <f> is given by 



ml{t) 



m. 



(24) 



The remarkable thing is that this is solvable by a Gaussian, even when the mass is time 
dependent. 



i 



-4> 2 /2a{tf 



27ra(t) 



(25) 



Here a it) obeys the equation 



Define S(t) = a 2 {t). Then 



da(t) 
dt 



dS(t) 
dt 



H 3 



in- 



'■a{t) 



(26) 



H3 + 2 -^s it) 



Air 2 



3H 



(27) 



Notice that this equation is readily interpreted as the field <f> subject to the force from the 
classical potential (the second term) in addition to the force driving Brownian motion due 
to de Sitter fluctuations [Ej]] (the first term). This equation is readily solved by finding 
the appropriate integrating factor and imposing the boundary condition that at t — — oo, 
a(— oo) = 0. The solution is 



H 2 r N 1 



-4/^ AT/3 



4/i 



*(N-N') 



dN' 



(28) 



where £ = fly/ [if and w = 1/ fi^fj,^. The integrand has a peak at N' = 0, with a width 
of order w. When £ is small, which is generally the case in our models, the peak is nearly 



14 



Gaussian and a saddle point approximation becomes applicable. So for N somewhat bigger 
than w (so that the peak is covered by the integration) and £ -C 1, S(N) is well-approximated 
by 

where the final approximation is valid if w < iV 

It should be borne in mind however that the Fokker Planck equation we used incorporated 
the slow roll equation of motion, which is valid for N small compared to w 2 Because it will 
turn out w is not large, the Fokker Planck description will be valid only at early times. We 
therefore use the Fokker Planck equation to establish initial conditions, and then use the 
classical field equations to describe the <fi evolution which we evolve numerically. 

While Eq. provides an analytic solution to the differential equation (E7p, the quali- 
tative behavior of the solution can be seen by looking at the differential equation itself. For 
large negative values of t, rn^it) is large and negative, providing a strong restoring force. 
This period is characterized by quasi-equilibrium evolution, in which the restoring force holds 
S(t) very close to its equilibrium value, —3H A /87r 2 m 2 i (t), for which dS/dt would vanish. The 
spread of this equilibrium probability distribution approaches zero in the asymptotic past, 
and grows monotonically with time. As t approaches 0, however, this equibrium value of 
S diverges, and the quasi-equilibrium regime ends because S(t) is not able to keep up. We 
can estimate when the quasi-equilibrium regime ends by asking when the velocity of the 
equilibrium value exceeds the diffusive velocity, H 3 /4ir 2 . For small £, this happens when 

N = Ni = — w-y/9/8. Then S(t) starts to grow diffusively, increasing linearly in time ac- 
cording to the first term on the right-hand side of Eq. (|27|). Neglecting the growth before 
N = Ni, which, in practice, changes the result by a number of order unity, we estimate 
S(0) as (H 2 /4:Tr 2 )(N — Ni) = \j9/8 H 2 w/An 2 , which in the limit of small £ gives an answer 

a factor of \/2ti smaller than the exact solution. 

The diffusive regime ends when the second term on the right-hand side of Eq. (p7|) be- 
comes larger than the diffusive term. This final phase can be called the classical regime, since 
the second term represents purely classical evolution. If only this term were included, the 
Fokker Planck equation would describe an ensemble of classical trajectories. This classical 
behavior is essential to our treatment of the problem, since it allows the description at late 
times to join smoothly to the full classical equations of motion which remain valid outside the 
slow-roll regime. The transition from the diffusive to the classical regime can be estimated by 
the "velocity matching criterion" , which is precisely when the two terms on the right-hand 
side of Eq. ( p7|) are equal, approximating the solution until this transition by the diffusive 
relation S(t) ~ H 3 t/4iT 2 . In the limit of small £, this velocity-matching condition holds at 

N = N Q = w^/9/8, and the value of the spread is given by S(N ) = <fi 2 = 3\/2H 2 w/16tt 2 . 
The classical regime can be approximated by constructing a solution to the classical equa- 
tions for 4>(t), starting from the initial condition <p(N ) = fa. If the asymptotic behavior of 
this classical solution (in slow-roll approximation) is compared with the asymptotic behav- 
ior of \J S(t) as given by Eq. (p9|) , it is found to be smaller by a factor of (Sne) 1 ^ ~ 2.9. 
In practice, we use the Fokker-Planck equation to establish the initial condition at iVo-In 
our numerical calculations we corrected for this discrepancy by using the initial condition 
(j)(N ) = & = (Svre) 1 /^. 
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We have determined the time evolution of (p and ip subsequent to the velocity matching 
time numerically. However, as for ip, the classical evolution of <p can be determined very 
well analytically. Again, we have to divide the analysis into three stages, according to the 
behavior of ip. 

At early times, the equation of motion for (p is approximately given by 

d* = ^4 (30) 

dN 9 y K ' 

which is solved by 

<j,(N) = fceW-Ntfl** (31) 

where we have imposed the boundary condition (p = (pi at N = N . This solution has 
assumed slow-roll which is only approximately valid. This stage of evolution of (p lasts for 



Ni ss 9/2y log(0 c /0j) e-folds, where (p c = \/lM'H. Numerically, we have found this 
answer to be off by 1-3 e-folds due to the correction to slow roll. 

At later time, as discussed above, ip begins to oscillate. Depending on the ip decay rate, 
there can be several e-folds between this time and the time at which the ip field decays. 
During this range of time, it can be checked that the ip 4 (p 2 term is no longer important to 
the equation of motion and that (p essentially keeps up with the minimum of the potential 

V = -=£ + *4 (32) 
2 8M' 2 v ; 

and is 

= = ^V^> /2 = /^7 e s^/ 2 (33) 

ip(N) ip c v 

which is valid when r <C H. Finally, the ip field decays. This occurs when T m H. For r 6 , 
the number of e-folds during this stage is approximately iV| = 0. For Ti, the value of the <p 
field when ip decays is approximately (pi = (H M' 3 'M£ j 'n) 1 / 6 . The total number of e-folds in 
this stage is approximately 2/3 log(0^/0 c ) which is N l 2 = (1/9) log(Mp/8if 2 7r). 
After ip decays, <p follows the ip = equation of motion according to 

(Pit) oc e rN (34) 

where 

/ n q 

(35) 

The number of e-folds in this stage is N 3 = (1/r) \og(<pf/{<p c , (pi}) depending on the decay 
rate, where (pf = ix/\/2f. 

In fact we have checked that the solution above gives the number of e-folds for inflation 
to end correct to within 1-3 e-folds. 

The reason we require an accurate determination of the number of e-foldings required 
for inflation to end is that it must be that the density fluctuations relevant for the observed 
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physical scales have exited the horizon while the field is in the early quasi-equilibrium stage. 
As we will now show, the quantum fluctuations of the field during the diffusive regime 
generate a spike in the density perturbation. For the model to be viable, it is important 
that this spike occurs at short wavelengths, so as to avoid conflict with observations. 

One interpretation of the source of density fluctuations is that the end of inflation occurs 
at different times at different points in space. The time delay function Ar(x), multiplied 
by H, is of the order of Sp/p at the time the wavelength re-enters the Hubble length. For 



example, if if) ~ e 



-m%N/3H 



and Aip ~ H/2tt, the density fluctuation constraint would be 



HAr = H 



6 x 10" 



Ai/j 1 



(36) 



Since has a more complicated evolution, the calculation of density fluctuations at early 
times is subtle. The problem is to estimate the density fluctuations caused by quantum 
fluctuations in ||21|| , whose value at times near t = is determined by the Fokker-Planck 
diffusion equation. Let us first consider the density fluctuations in early stages of the 
evolution. Suppose at a given time ti the solution P(0, t) is modified by displacing the 
entire probability distribution by an amount A0, so that P new (4>,ti) = P id(4 ) ~ A0, ii). We 
now ask how this will affect the time at which inflation ends. Because the if> and equations 
of motion are effectively decoupled, we can treat the if) field as uninfluenced, and treat the 
field as a free field evolving with time-dependent squared mass given by Eq. (p4j). 

We guess the solution is a shifted Gaussian, 



P(0,t) 



-{cf>-W/Mt? 



2na(t) 



(37) 



We find that the Fokker-Planck equation is satisfied provided that u 2 (t) obeys Eq. (|27j), and 
4>(t) obeys the equation 



dt 



3H 



(38) 



Imposing the initial condition 0(iVi) = A0, the solution to this equation is 



(f)(N) = A0exp 



4£ 2 



-4^iV/3 



-4^JVi/3 



+ -//J (N - N t ) 



(39) 



Since the entire distribution is shifted uniformly by 0(t), the implication is that so is each 
of the trajectories in the ensemble. The generic classical trajectory is the one whose value 
is equal to the RMS value of the distribution at large times as given by Eq. (28). Now by 
setting At = A0/0, we find 



HAt 



2^A0 /9tt 
Hptj, ' 



3\ 1/4 



exp 



(40) 



Here Nfi denotes the time that the initial condition is established, as the wave goes outside 
the Hubble length during inflation, and N e denotes the time at which inflation ends; both 
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times are measured in units of e-foldings. This formula applies for the first type of models; 
the exponent has a 2fi 2 b N/3 in the second type. 

The effect of a fluctuation is to determine the time and the value of at which diffusion 
ends and the classical evolution takes over (N e ). The direct change in the time at which 
classical evolution begins translates into a difference in time at which inflation ends. We 
expand the above answer for small £ to get 

HAr = (9nY 4 ——w 3 / 2 e- 2N2 f^ 2 (41) 
2N e H 

where Nfi is the time at which the fluctuation occurs. This falls off from the peak at Nfi = 



like a Gaussian with width J 9w 2 /2. What this tells us is that fluctuations formed sufficiently 
early (or late) will not delay the onset of the classical regime significantly and not give a 
significant contribution to density fluctuations. However, fluctuations formed during the 
diffusive growth regime near N = will create far too large density perturbations. These 
fluctuations must be such that they are not relevant to observable scales. We can observe 
back to about 40 e-folds before the end of inflation, so we require that the fluctuations 
formed at this time were sufficiently small, or that inflation must end by 40 e-folds beyond 
the time when the fluctuations satisfy the experimental bound. From Eqn. II, one can 
deduce this time is approximately N d f = — 5u>. It might be thought that another solution 
is that inflation ends very slowly, so that 40 e-folds before the end one is in the classical 
regime. However, when is evolving according to the classical equations of motion, the 
scale dependence of density perturbations is much too large. 

So the number of e-folds beyond the time when density fluctuations in are sufficiently 
small is 

[g 

N=\N df \ + \-w + N 1 + N 2 + N 3 (42) 
V 8 

Density fluctuations on the scale of 1 Mpc are formed 

N Mpc = 38 + | log(M/10 n GeV) + ^ log(7W/10 7 GeV) (43) 

e-folds before the end of inflation. We require that iV is less than Nm pc - By following 
through the above calculations, one can see this gives the approximate constraint w ~ 1. 
The detailed application of the constraint gives the constraints illustrated in Figures 1-3. 
These plots were made assuming the larger decay rate. The total number of e-folds for 
the same parameters is generally about 5 larger with Ti which can be accomodated with a 
modest change in //^. 



VII. BLACK HOLES? 

Because of the large peak in the density perturbation spectrum on small length scales 
arising from the contribution, there is a danger too many small black holes being created. 
There are fairly strong constraints on the fractional mass density in black holes on small 
scales Hl5| . We investigate these constraints on our model in this section. 
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First we summarize the constraints. In the paper of Carr, Gilbert, and Lidsey, the 
constraints are presented in several forms; one constraint is on the parameter 5 which is 
related to 5p/ p. In the mass range above 10 30 gm there are bounds from CMB distortions 
constraining 5 to be less than about 10~ 2 . In the mass range between 10 30 gm and 10 10 gm 



the bound on 5 is approximately 10 . In |15[ a bound due to relics is deduced constraining 



S between about 10 -2 and 10 -L5 in this mass range. This constraint from relics is perhaps 
more speculative than robust bounds from not exceeding critical density, or that decaying 
black holes do not produce too much entropy. 

In order to apply these bounds, one needs to know the probability of black hole creation 



as a function of 5. Based on [22], the bound is obtained from applying the formula for the 



probability of a region of mass M forming a primordial black hole 

ft(M)«<f(M)exp(-i) (44) 

where the equation of state when the perturbation enters the horizon is p = jp. For the 
scales which are of interest to us, 7 = 1/3, which we will assume in the equations below. 

The above bound comes from considering a spherically symmetric overdense region. The 
requirement is made that when the overdense region stops expanding, it's size S c exceeds 
the Jeans radius Rj at this time t c , in order to collapse against the pressure. To derive 
exact numerical bounds on 5 requires that this be the precise condition. Without solving 
the full problem explicitly including the pressure effects near the boundary it is difficult to 
state precisely the conclusion, which gives rise to some overall uncertainty in the 5 bound. 
However, one should be able to obtain a conservative bound on 5 through this approximation. 
However, even using this approximation, we find numerical discrepancies with the precise 
production rate which would be predicted. First, the relation between t c and £ (the time 
at which the perturbation begins to evolve separately from the homogeneous background 
in which it is embedded) should be t c = 2t 5~ 1 (the 2 being omitted in Ref. p2| ) and the 



relation between Rj and t c should be Rj = 2v s iT^2/3t c rather than v s 2t c (in actual fact the 
2 was omitted but cancels later on), where v s is the sound velocity. Overall this translates 

into the bound 5o > v 2 2n^j2/3 (jf^ (where the correct relation S c = Ro8 1//2 has been 
substituted and the ratio t /R has been replaced by the appropriate mass ratio). The 
implication is that the factor 7 2 /2 in Eqn. [0] should be replaced by 47r 2 7 2 /3, which in turn 



decreases the strength of the bound on 5 by a factor of about 5. This will of course also 



weaken the bound on the scalar index n given in [T3[ . 

Our spectrum is not scale invariant on these small length scales which is important when 
calculating S from HAt. However a very conservative upper bound on our spectrum is a 
scale invariant spectrum starting at a small length scale (near the peak of the Gaussian) 
and which is constant over smaller wavelengths. This spectrum would be a scale invariant 
spectrum with a cutoff at large wavelength, and can readily be compared with the analysis 



of Ref. | 15fl . The normalization they use for 5 can be extracted from their Eqs. (4.2)-(4.4), 
which express the value of 5 at the COBE scale in terms of the underlying inflaton potential. 
Assuming that T/S <C 1, their equations reduce to 5 = 0.99V 3 / 2 /V . By comparing with 
the relation HAt = H5<j)/4> = H 2 / (27T0) , one finds that S can be related to the fluctuations 
in r by 5 = .086HAt. We can then apply formula ^0] to find that 5 (normalized as above) 
never exceeds 0.02, in the parameter regime presented in Figures 1-4. Because the quoted 
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constraint on 5 in [15| appears be too strong by a factor of 5, we conclude that we are 
consistent with reasonable estimates of the black hole constraint. Therefore, even with a 
conservative overestimate of S, the constraint from black holes is satisfied. However, there 
can be a sizable fraction of matter in black holes, which would be interesting to study in the 
future. 

As a final comment, we remark that the bound from black holes is somewhat weaker if 
inflation ends more quickly, which it does for lower w. In fact, for w ~ 0.5 (corresponding to 

about four times larger than assumed) inflation ends in about 10 e-foldings. In this case, 
only the bounds from relics would apply. As this bound is more speculative, it is possible 
that even large perturbations on this scale would be acceptable. In reality, smaller w, while 
decreasing the length of time for inflation to end, also decreases the density perturbations. 
Larger always leads to a smaller fraction of the universe in black holes. 

In summary, the black hole constraint is a serious constraint and must be accounted for. 
This is another constraint which would forbid large w, since the maximum value of HAt 
grows as w 1 ^ 2 (when the dependence of N e on w is accounted for). However we have seen 
our model is safely within the bounds given in the literature once we have applied the 
bound given in Section 7. 

However these bounds are not sufficiently precise at present and it would be interesting 
to do a more accurate calculation of the mass fraction in black holes both for our model and 
in general.^ The effects we have discussed should weaken existing bounds. 

VIII. GRAVITINO CONSTRAINT 

In this two field model of inflation, the source of entropy and energy in the universe is 
the decay of the ip an d fields. The ip field decays first, as discussed earlier. We assume 
the decay products are quickly thermalized, giving an effective temperature T^. Sometime 
afterward the field reaches the minimum of its potential and begins to oscillate about it. We 
assume that these oscillations are damped by gauge or Yukawa couplings. As in Ref . , the 
decay can occur through a coupling g 2 / (0) / d 4 0x [ X ( t ) or through a direct Yukawa coupling 
A00i0 2 . This leads to a reheat temperature equal to max(g 2 ^ 3 m 5 J 6 M 1 / 6 , m^/VA"), which is 
generally of order 10 5 -10 7 GeV. Since most of the energy of the universe evolves from the 
coherent oscillations of the <fi field, this reheat temperature sets the initial conditions for the 
subsequent evolution. As discussed in Ref. p3| , this reheat temperature is low enough to 



avoid the overproduction of gravitinos, even if gravitinos are as light as 100 GeV. 

However, the initial temperature of the thermal plasma of ip decay products can be as 
high as 10 11 GeV, so the production of gravitinos by this plasma must be examined. In this 
section we show that this constraint is never more restrictive than the constraints already 
discussed. 

Gravitinos are produced by scattering processes of the thermal radiation, but interact at 
a rate suppressed by m^ 2 /M 2 . They are potentially dangerous since they are not thermal- 
ized and have a long lifetime. The most stringent bounds are obtained by considering the 



^We thank B. Carr for informing us that work is in progress on this subject. 
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influence of these late decays on nucleosynthesis. The exact bound depends on the gravitino 
mass, but for 100 GeV < m 3/2 < 1 TeV it is T R < 10 7 ~ 9 GeV ||. 



Neglecting decay when considering gravitino production, one writes the Boltzmann equa- 



tion for the gravitino number density as | 23fl 

+ 3Hn 3/2 = a cS n 2 ad , (45) 

where n TSu d is the equilibrium number density of a single species of scalar boson (n rad = 
(C(3)/7r 2 )T 3 ), and both the cross sections and the multiplicity of species are accounted for 
by the factor a e s. In terms of Y 3 / 2 = n 3 / 2 /n ra _d, we have 

—f- = a cS n rad . (46) 

To a reasonable approximation a c s can be taken as constant, although it does vary as the 
coupling constants run and as species freeze out from the thermal equilibrium mix p3j . For 
the standard case of a radiation-dominated universe, n rad oc 1/R 3 (t) oc 1/t 3 / 2 (where R(t) is 
the scale factor), so the total gravitino production can be estimated by integrating Eq. fl46|) 
from the initial reheat time to to infinity. This gives Y 3 / 2 ~ 2cr c g- n radi o to = cr c s ^rad,o/-^o 
To, where the subscript refers to the time of reheating. To obtain a reasonable estimate of 
the present value of Y 3 / 2 , one must divide this value by a dilution factor to account for the 



production of photons at times much later than to- According to Ref. [23], the final result is 



Y 3/2 (T « IMeV) « 2.14 x 10" 11 ( ^J^ ) • (47) 

To derive a conservative estimate for the gravitino production of the if) decay products, 
we assume that the energy released by the ip decay is approximately equal to the energy 
stored in the oscillating field when inflation ends. According to the numerical simulations 
of our model, the fraction of energy in the if) field was generally less than this by a factor 
of at least 10 4 , except in the case M' = M v and a slow decay rate, in which case ip can 
store a substantial fraction of the energy at the end of inflation. The universe then rapidly 
become matter-dominated, so R(t) oc t 2 / 3 . Repeating the calculation for Y 3 / 2 with this time 
evolution, one finds Y 3 / 2 = | a e fj n Tad /H , essentially the same formula as above. 

However, in this model there is an additional dilution of the ip decay products, because 
the field behaves as a coherent state of nonrelativistic particles for a time T^" 1 , and then 
the 4> particles decay to produce radiation. Before the particles decay, the energy density 
Pip of the ip decay particles (assumed to be effectively massless) is suppressed relative to the 
energy density of the field by one power of the growth of the scale factor between the 
times to and T^ 1 , which is l/(r^to) 2//3 . When the particles decay to radiation, the number 
of radiation particles produced exceeds the number of ip decay particles by (p^/p,/,) 3 / 4 ~ 
l/^T^to) 1 / 2 . Relating to the reheat temperature Tr ~ 10 7 GeV and taking t ~ 
the dilution factor is found to be approximately (Mpm^/T^) 1 ^ 2 w 10 4 . Incorporating this 
extra dilution factor into Eq. (f47l), we find that gravitino production from if) decay products 
give 

Y 3/2 (T « 1 MeV) - 10- 15 L^^) ■ (48) 
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where we have explicitly incorporated our assumption that ip and initially carry comparable 
energy. 

Thus, a if) reheat temperature of 10 11 GeV produces no more gravitinos then a final reheat 
temperature of 10 7 GeV. Thus, we find that no further constraints need to be imposed. 



IX. BARYOGENESIS 

In the context of late scale inflation, it is worthwhile to investigate the question of how 
baryons are created. There are essentially two known possibilities. Electroweak baryoge- 
nesis |24| is possible since the reheat temperature will generally be above the weak scale. 
Alternatively, a model in the context of supersymmetry invites investigation into the Affleck 
Dine scenario p5 |. 



In the Affleck Dine baryogenesis scheme, a field x which carries baryon number acquires 
a large displacement relative to its true minimum somewhere during the early evolution of 
the universe. If the interactions which drive the field to the true minimum are CP and 
B violating, the field will store baryon number, and subsequently decay to baryon number 
carrying particles. 

In our model, in principle, the fields <fi or if) could be the AD fields. However this does 
not work. The problem is that fields which carry baryon number will generally also carry 
charge, so that <f> is not a good possibility since charge (or color) would be spontaneously 
broken by the vacuum. Although if) is in principle a candidate, the ratio of baryon number 
stored by the if) field to entropy will be too small. 

This can be deduced from a detailed study of the if) field. The first point to observe 
is that the potentials we have studied to now are B and CP conserving. This is because 
we have neglected the soft U A" type terms and possible cross terms which can violate CP. 
When these are included, we find there can cause a small change in the detailed evolution of 
the field. At the time the A (CP and B violating terms) are large, the if) field only carried 
a small fraction of the energy of the universe. 

However, a separate flat direction which plays the role of the AD field would work. If 
the AD field is independent of the inflation fields, it can then have large expectation value 
through the final stages of inflation. If H is somewhat larger than itlad, the analysis is 
similar to that in Ref. [[16] where a much larger H was assumed. The baryon to entropy 
ratio is approximately 

s n AD m AD T£ 

where n^/nAD gives the baryon to particle number ratio in the AD field, and should be 
order unity if the potential for the field is B and CP violating. The last factor pad/Tr is 
determined by the amplitude of the AD field at the time it evolves towards its true minimum, 



which is determined by higher dimension operators in the potential [16]. One can readily 
obtain acceptable values for the baryon density if the dimension of the operator in the 
superpotential which lifts the flat direction is greater than 4. The lower reheat temperature 
expected in these models requires a correspondingly larger factor pad/T^, so a dimension 4 
operator in the superpotential which lifts the AD field is insufficient. 
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X. DISCUSSION 



There are several comments to make about the models we have considered. First there 
is the fact that we are considering very late inflation. We do not address the question of 
why the universe has lasted to this point |J. We have only addressed the question of a 
late inflationary epoch responsible for solving the horizon and flatness problems and for 
generating the necessary density fluctuations. 

One possible solution is that the initial ijj value exceeds M p . In this case, it is possible 
that chaotic inflation could solve the problem raised above. However, subsequent to this 
stage of inflation, one would expect inflation as described in this paper which would create 
density fluctuations of the right size. 

Another point we have not addressed is what our models look like when embedded in 
supergravity theories. Our point of view throughout this paper is to regard the theory as 
an effective theory expanded in powers of M p (and M'). We have neglected terms which 
are suppressed by higher powers of the Planck scale. For the same reason we have assumed 
a minimal Kahler potential when deriving the potential. From this point of view, any 
secondary minima which occur in supergravity for field values exceeding M p are not to be 
trusted. 

One aspect of our models which is important is the requirement of renormalizable cou- 
plings of the fields in the direction in order to obtain sufficiently high reheat temperature. 
For this reason we expect it is more likely that the fields and ip correspond to flat directions 
of a renormalizable theory (along the lines discussed in Section 6) than to true moduli fields 
(of string theory) [^(| . 

In our models, we saw that there was usually some small but not very small number. 
Either M' is M p in which case ^ and 1//^ are of order 100, or M' is smaller than M p which 
permits and fi^ closer to unity. Another possibility is that the small number is related 
to a Yukawa coupling. We discuss each of these possibilities in turn. 

It is well known that there can be H dependent correction to the soft supersymmetry 
breaking masses at early times when H exceeds m 3 / 2 [|l6|,^7j]. This means that small fi^p is 
necessarily obtained by tuning. Stewart [123 has presented criteria which are sufficient for 
the cancellation of supergravity corrections to the inflaton mass, so that even a large ratio is 
technically consistent. However these conditions will only work when the scale of inflation is 
above the supersymmetry breaking scale. One therefore needs to invoke a new mass scale. 
Generally Stewart chooses the scale of gaugino condensation. It is hard to see how this scale 
is realized in an actual model although it could present an interesting alternative. There is a 
tradeoff between the complexity of the model and the "naturalness" of taking somewhat 
smaller than H. 

In the models where fi^ is small, we found that consistency of the model required that 
Hj> is large, with the product [i<$,\x^ being a number of order unity. It might be thought that 
this large value of /i^, could be explained as due to large H dependent mass corrections. 
However, it is not possible to introduce a large /i^ without the tuning parameter appearing 
in some other unnatural feature of the potential. For example, a large H dependent mass 
could introduce a new minimum for <ft which is closer to the origin so that the VEV of <fi is 
correspondingly smaller than M p . 

We have seen however that the tuning of mass ratios is significantly reduced if we accept 



23 



a higher scale than the conventional intermediate scale as setting the overall energy density. 
However if this scale has anything to do with visible supersymmetry breaking, the highest 
scale possible is probably the gaugino condensate scale. The predictions for n and R would 
be very similar, so the general test for this class of models would still be valid. 

We regard the small tuning of parameters as a necessary aspect of the models with 
M' ~ M p . The necessary numbers may or may not be present. The tuning is certainly 
much smaller than in a typical inflation model. 

On the other hand, M' might be smaller. This requires the presence of another mass 
scale in the theory. If this lower scale exists, one can obtain fi^ and fi^ closer to unity. 

The other possibility is that there is a small Yukawa coupling. This is probably not 
such a bad possibility. First of all, the necessary coupling is no smaller than known Yukawa 
couplings and might even be related to them as in the model of Section 6 and obvious 
variants. Second, known Yukawas can be derived as the ratio of mass scales. In the SU(6) 
model we discussed this has been done in Ref . |2T||H| . It is not unreasonable to think there 
might be such effective Yukawa couplings in hidden sectors of the theory, as well as in the 
single known sector. 

In most known hybrid inflation models other than the one we discussed, the ip field is 
very light, while the field is very heavy, of order M. This is a more serious technical 
problem since radiative corrections will generally give ip too large a mass 0. Even if the 
model is supersymmetric, supersymmetry breaking during inflation would induce a large 
mass for tp. Although at tree level Str(M 2 ) = 0, this is not sufficient to prevent radiative 
corrections at higher loop order. One can perhaps allow for such a hierarchy, but at the 
expense of additional complexity and mass scales. A chief advantage of our model is that 
both and are of order the soft supersymmetry breaking scale so this problem does 
not arise. 

We view our model as the simplest illustration that flat directions of supersymmetric 
theories are consistent with the requirements of inflation when one allows for more than one 
field in the inflation sector. It is likely that the small parameters which might be required 
(of order 0.01 to 0.1) are present. Alternatively there might be more subtle mechanisms at 
work. Either way, one would conclude that the scale of inflation is very low. Even allowing 
inflation to be determined by the higher gaugino condensate scale, one would conclude that 
H during inflation is between 10 4 and 10 7 GeV, and tensor perturbations are small. 

It is important that there are observational consequences to this type of model. The 
combination of measuring the scalar index n and T/S should either rule out or encourage 
belief in the mechanism at work here. As discussed in Section 3, 



Because the second term is negligible in models of the sort we are considering, where i/j at 
the end of inflation is much less than Mpi, it is only the last term which causes the deviation 
of n from unity. If the dependence on ip is dominated by a mass term, as in the model of 
Section 3, the correction to n will be positive (but small). We then expect n greater than 
or equal to unity, and T / S to be small. 




(50) 
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XI. CONCLUSION 

We have shown that with more than one field it is possible to construct models of inflation 
with no small parameters. Furthermore, the mass scales which seem to most naturally 
appear in these models are of order m 3 / 2 , about 1 TeV, and Mj, about 10 11 GeV, leading to 
a natural association with supersymmetric models. These models give rise to the correctly 
normalized density perturbations, even though the Hubble constant is quite low, of order 
10 4 GeV, because the value of the inflaton field at the end of inflation is much lower than 
the Planck scale. The key to producing more such models is a sensitive dependence of the 
4> potential on the value of the ip field, so that the motion of the ip field can trigger the end 
of inflation while its value is small. 

It seems that multifield models are probably the most natural models which can imple- 
ment inflation with weak scale Hubble constant, and that furthermore, these are probably 
the most natural inflation models in that they involve no new small parameters. The requi- 
site small parameters arise naturally from the ratio of mass scales. These models have the 
further advantages that they can be explicitly realized and one can calculate the relevant 
parameters for any particular implementation. They might even occur in simple extensions 
of the MSSM. 

Perhaps the most important property of a model is its testability, and our proposed 
models have several characteristics that are in principle observable. The scalar index n 
which characterizes the scale dependence of density perturbations is always greater than 
unity. It is very close to unity for the model of Sec. Ill with M' at the Planck or GUT scale, 
but for M' at the intermediate scale or for the model of Sec. IV, it could be as large as 1.2 
for the parameters shown in our plots. In all cases tensor perturbations are negligible. An 
especially distinctive feature is a large spike in the density perturbation spectrum at present 
wavelengths of about 1 Mpc or less. 

Acknowledgements 

We are very grateful to Andrew Liddle, David Lyth, and Ewan Stewart for discussions 
and the Aspen Center for Physics where these discussions took place. We also thank Sean 
Carroll, Csaba Csaki, Arthur Kosowsky, Andrei Linde, and Bharat Ratra for their comments. 
We thank Bernard Carr, Jim Lidsey, Avi Loeb, Paul Schechter, and Paul Steinhardt for 
discussions and correspondence about black holes. We also thank Krishna Rajagopal for his 
comments on the manuscript. 



25 



REFERENCES 



[1] For a review, see A. D. Linde, Particle physics and inflationary cosmology (Harwood 

Academic, Switzerland, 1990). 
[2] A. Albrecht, S. Dimopoulos, W. Fischler, E.W. Kolb, S. Raby, & P.J. Steinhardt, Nucl. 

Phys. B229, 528 (1983) 
[3] K. Freese, J. Frieman, and A. Olinto, Phys. Rev. Lett. 65 (1990) 3233. 
[4] T. Banks, M. Berkooz, P. Steinhardt, Phys. Rev. D 52 (1995) 75. 
[5] T. Banks, M. Berkooz, G. Moore, S. Shenker, and P. Steinhardt, preprint RU-94-93, 

[hep-th 9503114| . 
[6] S. Thomas, Phys. Lett. B. 351 (1995) 424. 

[7] G. Dvali, Q. Shafi, R. Schaefer , Phys. Rev. Lett 73 (1994) 1886. 
[8] L. Randall, C. Csaki, PASCOS Proceedings (1995), |hep-ph 9508208 



[9] L. Randall and S. Thomas, Nucl. Phys. B 449, 229 (1995). 
[10] A. D. Linde, Phys. Lett. B259 (1991) 38; A. R. Liddle and D. H. Lyth, Phys. Rep. 231 

(1993) 1; A. D. Linde, Phys. Rev. D49 (1994) 748; E. J. Copeland, A. R. Liddle, D. 

H. Lyth, E. D. Stewart and D. Wands, Phys. Rev. D49 (1994) 6410; E. Stewart, Phys. 

Lett. B 345,414 (1995). 
[11] A. Linde, Phys. Lett. B 129, 177 (1983). 

[12] A. Guth and S.-Y. Pi, Phys. Rev. Lett. 49, 1110 (1982); A. Starobinsky, Phys. Lett. B 

117, 175 (1982); S. Hawking, Phys. Lett. B 115, 295 (1982); J. Bardeen, P. Steinhardt, 

and M. Turner, Phys. Rev. D 28, 679 (1983). 
[13] A. R. Liddle and D. H. Lyth, |astro-ph/ 9409077 . 
[14] M. Turner, Phys. Rev. D48 (1993) 5539. 
[15] B. Carr, J. Gilbert, J. Lidsey, Phys. Rev. D 50 (1994) 4853. 
[16] M. Dine, L. Randall, and S. Thomas, Phys. Rev. Lett. 75, 398 (1995). 
[17] T. Gherghetta, C. Kolda , S. Martin, |hep-ph/9510370| (1995).. 
[18] K. Inoue, A. Kakuto and H. Takano, Prog. Theor. Phys. 75 (1986), 664; 

A. Anselm and A. Johansen, Phys. Lett. B200 (1988), 331; 

A. Anselm, Sov. Phys. JETP 67 (1988), 663; 

R. Barbieri, G. Dvali, A. Strumia, Nucl. Phys. B391 (1993), 487; 

Z. Berezhiani and G. Dvali, Sov. Phys. Lebedev Inst. Rep. 5 (1989), 55; 

R. Barbieri, G. Dvali, M. Moretti, Phys. Lett. B312 (1993), 137. 
[19] Z. Berezhiani , C. Csaki, L. Randall., Nucl. Phys. B 444 (1995) 444. 
[20] A. Linde and A. Mezhulmian, |gr-qc/93040T5] 

[21] T. Bunch and P. Davies, Proc. R. Soc. A 360, 117 (1978); A. Linde, Phys. Lett. B 116, 
335 (1982); A. Starobinsky, Phys. Lett. B 117, 175 (1982); A. Linde, Phys. Lett. B 131, 
330 (1983). 

[22] B. Carr, Ap. J. 201 (1975) 1. 

[23] T. Moroi, Ph.D. Thesis |astro-ph/9503210| , and refs. therein. 

[24] For a review see A. Cohen, D. Kaplan, and A. Nelson, Ann. Rev. Part. Nucl. Sci., 43, 
27 (1993). 

[25] I. Affleck and M. Dine, Nucl. Phys. B 249, 361 (1985). 

[26] L. Dixon, in Superstrings, Unified Theories, and Cosmology 1987, Proceedings of the 
Summer Workshop in High energy Physics and Cosmology, eds. G. Furlan et. al, (World 
Scientific, Singapore, 1987). 



26 



[27] M. Dine, W. Fischler, D. Nemeschansky, Phys. Lett. B 136 (1984) 169; 

O. Bertolami and G. Ross, Phys. Lett. B183 (1987)163; 

E. Copeland et. al. Ref. Jig. 
[28] E. Stewart, Phys. Rev. D 51 (1995) 6847. 

[29] R. Barbieri, G. Dvali, A. Strumia, Z. Berezhiani, L. Hall, Nucl. Phys. B432 (1994), 49. 



27 



